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ABSTRACT 


W a  developed  Micro  SAINT  computational  networks  for  numerical  integra¬ 
tion  and  solving  initial  value  problems  for  linear  and  nonlinear  first-  and 
second-order  ordinary  differential  equations  as  well  as  for  systems  of 
differential  equations.  These  Micro  SAINT  computer  programs  are  written 
with  a  user  friendly  approach  where  the  user  will  be  required  to  supply  the 
input  information  and  the  functional  form(s)  of  the  function(s)  in  the 
"function  library"  section  of  Micro  SAINT  without  any  changes  i:i  the  main 
programs. 

These  computational  modules  could  be  used  as  subnetworks  in  modeling 
psychophysiological  and  biomedical  problems  of  interest  in  naval  aerospace 
medical  research.  For  example,  Micro  SAINT  developed  models  can  be  used  by 
staff  medical  officers  to  predict  psychophysiological  performance  of  naval 
aircrew  personnel  under  sustained  operational  work  schedules. 
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1,  INTRODUCTION 


In  a  2-day  workshop  on  Micro  SAINT  analysis  given  by  S.  E.  Shamma  on 
November  10  and  17,  1988,  for  staff  members  of  the  Naval  Aerospace  Medical 
Research  Laboratory,  Pensacola,  Florida,  a  question  was  raised  about  .the 
potential  of  using  Micro  SAINT  software  [1]  for  numerical  solutions  of 
mathematical  problems,  especially  differential  equations.  In  this  report, 
we  answer  this  question  affirmatively. 

We  developed  Micro  SAINT  networks  for  numerical  integration  and  solv¬ 
ing  initial-value  problems  for  linear  and  nonlinear  first-  and  second-order 
ordinary  differential  equations  as  well  as  for  systems  of  differential 
equations.  These  computational  modules  are  expected  to  be  used  as  subnet¬ 
works  in  modeling  problems  of  interest  in  naval  aerospace  medical  research. 

The  following  sections  constitute  a  detailed  theoretical  summary  and 
Micro  SAINT  programs  for: 

a.  Numerical  integration  of  a  function  f(x)  using  Trapezoidal  Rule. 

b.  Numerical  integration  of  a  function  f(x)  using  Simpson's  Rule, 

c.  Composite  Simpson's  Rule  for  double  integrals. 

d.  Composite  Simpson's  Rule  for  triple  integrals. 

e.  Euler  method  for  solving  a  first-order  ordinary  differential 
equation. 

f.  Modified  Euler  method  for  solving  a  first-order  ordinary  differen¬ 
tial  equation. 

g.  Runga-Kutta  method  of  order  four  for  solving  a  first-order 
ordinary  differential  equation. 

h.  Runga-Kutta  method  of  order  four  for  solving  a  first-order  system 
of  ordinary  differential  equations. 

i.  Runga-Kutta  method  of  order  four  for  solving  second-order  (linear 
or  nonlinear)  ordinary  differential  equations. 
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The  Micro  SAINT  computer  program;  are  written  with  a  user  friendly 
approach.  The  user  will  supply  the  input  information  in  specified  places 
in  the  program,  and  the  output  is  stored  in  the  snapshots  output  files.  The 
input  information  shall  be  entered  in  the  "function  library"  section  of 
Micro  SAINT.  The  initial  and  end  conditions  shall  be  entered  in  a  function 
named  "initl"  (initial).  It  consists  of: 

intilx  =  initial  value  of  the  independent  variable  x, 
endx  =  end  value  of  the  independent  variable  x, 

numir.cvl  =  number  of  subintervals  for  integration  or  subdivisions  in 
the  case  of  differential  equations;  and  initial  y  or  initial  u^  and  U2  in 
the  case  of  solving  a  system  of  differential  equations. 

The  integration  function  or  the  functional  form(s)  of  the  differential 
equation(s)  shall  be  entered  in  a  straightforward  way  in  the  "function 
library"  of  Micro  SAINT.  We  will  illustrate  these  procedures  by  applying 
the  programs  on  examples  from  reference  [2].  The  main  computer  programs 
are  listed  in  appendix  A.  The  user  may  view  each  main  program  as  a  "black 
box"  since  the  input  information  is  entered  separately  in  the  "function 
library"  of  Micro  SAINT.  Some  results  are  presented  graphically  in  appen- 
d  ix  B. 

2.  SUMMARY  OF  NUMERICAL  INTEGRATION  METHODS  AMD  INPUT  INFORMATION 
TO  THE  "FUNCTION  LIBRARY*'  OF  MICRO  SAINT 

NUMERICAL  INTEGRATION 

We  considered  two  methods,  Trapezoidal  and  Simpson's  Rules,  for  comput- 
b 

f 

ing  the  integral J  a  f(x)  dx  and  a  composite  Simpson’s  Rule  for  double 
integrals. 
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1.  Trapezoidal  Rule 

If  £€  C2  [a,b]  with  h  ■  (b-a)/n  nnd  Xj  *■  a  jh  for  each 
,j  ■  0,l,2,„.,,n,  the  trapezoidal  rule  for  n  subintervals  is: 

b  -i 

J  a  f(x)  dx  -  _h[  ^  { f ( xj )  +  £(xj+1)}  ]  +  error, 

j-0 

Error  ■  0(h2). 

The  usor  shall  enter  the  input  information  in  the  "function 
library"  of  Micro  SAINT  as  follows: 

initix  **  a ,  endx  ■  b,  numintvl  »  n,  and  the  expression  for  the 
function  f.  The  following  example  illustrates  the  case  where  a  “  0,  b  ■  1, 
n  »  10,  and  f(x)  ■*  1  +  x. 

TA  . ,  Function  Library  Input  for  Intgtrpz  Program. 

FUNCTION  LIBRARY  Model  Name:  intgtrpz 

Name:  Expression: 

f  1+x; 

ini tl  numintvl-10; initlx-0;endx«l; 


2.  Simpson's  Rule 

If  f  €  C^[a,b],  with  h  -  (b»a)/n,  where  n  -  2m,  n  must  be  an  even 
integer,  with  x j  *  a  +  j  h  for  each  j  *  0, 1,2,,,. ,2m,  the  Simpson's  Rule  for 
n  subintervals  Is: 


f(x)  dx 


m 


m  £if(. 

I 


2  j-2 


_2^  C x2 j  —  1  ^ 


£<*2j)l} 


■f  Error , 


Error  -  0(h4). 

The  use::  shall  enter  the  input  information  in  the  "function 
library"  of  Micro  SAINT  aa  illustrated  in  the  trapezoidal  case. 
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3.  Composite  Simpson's  Rule  for  Double  Integrals 

The  program  "dblinteg"  approximates  the  double  integral  of  a  function 
f(x,y)  with  limits  of  Integration  from  a  to  b  for  x  and  from  c(x)  to  d(x) 
for  y,  using  a  composite  Simpson's  Rule.  To  evaluate  the  integral 
x  -  t 


/ 


y  -  d(x) 

r  f(x»y)  dydx, 


x  »  a  y  -  c(x) 

the  user  needs  to  supply  the  functions  f(x,y),  c(x),  and  d(x)  as  well  as  the 
parameters  initlx  =  a,  endx  «*  b,  and  the  number  of  divisions,  numdivx,  and 
numdivy,  along  the  x  and  y  axes;  "numdivx"  and  "numdivy"  must  be  even 
numbers.  The  approximate  value  of  the  integral  is  stored  in  the  output 
snapshot. 

As  an  illustration,  consider  the  integral 
x  “  1  y  -  x 

/  /  .  (xy)  dydx, 

x  -  0  y  -  x2 

Here  initlx  *  0,  endx  ■  1,  c(x)  *■  x  ,  d(x)  *»  x,  and  f(x,y)  «  xy.  The  user 
supplies  this  information  as  well  as  numdivx  and  numdivy  in  the  "function 
library"  of  Micro  SAINT  as  follows: 

TABLE  2.  Function  Library  Input  for  Dblinteg  Program. 


FUNCTION  LIBRARY 


Model  Name:  dblinteg 


Name : 
f 

initl 

funedofx 

funccofx 


Expression: 

x*y; 

ini  tlx*0;  endx*  1;  numdivx**  10;  numdivy*<10; 

x; 

x*x; 


The  exact  value  of  the  integral  is  1/24,  and  the  computed  answer  la  0.041650. 
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4 .  Composite  Simpson's  Rule  for  Triple  Integrals 


The  program  "triplint"  approximates  the  triple  integral  of  a 


and 


x  ■  b 

/  J 


limits  from  a  to 

b  for  x,  from  c(x)  to  d(x)  for  y. 

ft  (x,yv  for  z. 

To  evaluate  the  integral 

d(x)  z  » 

<x.y) 

/ 

f(x,y,z)  dzdydx, 

c(x)  Z  “ 

u,y) 

a  ■  y 

the  user  needs  to  supply  the  functions  f(x,y,z),  c(x) ,  d(x),  0(x,y), 

^(xjy)  as  well  as  the  parameters  ini tlx  -  a,  endx  *  b,  and  the  number  of 
divisions,  numdlvx,  numdlvy,  and  numdlvz,  along  x,  y,  and  z  axes, 
respectively;  "numdlvx,"  "numdlvy,"  and  "numdlvz"  must  be  even  numbers. 

The  approximate  value  of  the  Integral  is  stored  In  the  output  snapshot. 

As  an  illustration,  consider  the  integral 

x“l  y  ■  x  z  *  2 

f  J  J  (xyz)  dzdydx. 

x*  0  y-x  z  «  xy 

2 

Here  lnltlx  *  0,  endx  *  1,  c(x)  -  x  ,  d(xv  ■  x,  Qt  (x,y)  *  xy,  (x,y)  -  2, 
and  f(x,y,z)  *  xyz.  The  user  supplies  this  information  as  well  as  numdlvx, 
uuradivy,  and  numdlvz  in  the  "function  library"  of  Micro  SAINT  as  follows: 
TABLE  3.  Function  Library  Input  for  Triplint  Program. 


FUNCTION  LIBRARY  Model  Name:  triplint 


Name 

Expression: 

cx 

x*x; 

dx 

x; 

betaxy 

2; 

alphaxy 

x*y; 

fxyz 

x*y*z ; 

lnltl 

inltlx-0;endx»l;numdivx»10;riuradivy-10;numdlvz»105 

The  exact  value  of  the  integral  Is  0.078125,  and  the  computed  answer  is 
0.078578. 
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3.  NUMERICAL  SOLUTION  OF  DIFFERENTIAL  EQUATIONS 


a.  First-order  ordinary  differential  equation:  We  consider  three 
methods,  Euler,  modified  Euler,  and  Ruuga-Kutta  of  order  four  for  solving 
the  initial  value  problem. 

-  f (x,y) ,  y(x0)  -  y0. 

1.  Euler's  method: 

The  difference  equation  associated  with  Euler  method  is: 

yi  *  ^1-1  +  hf<xi-l*yi-l>  +  Error»  for  1  “  l,2,3,...,n,  where  n  -  number  of 
intervals,  x^  -  xQ  +  ih,  and  Error  ■  0(h2X 

The  user  shall  enter  the  input  information,  initlx,  endx,  numoivx, 
initly,  and  the  functional  form  of  f(x,y)  in  the  "function  library"  of 
Micro  SAINT. 

The  following  example  illustrates  the  input  information  needed  for 
solving 

dy  -  f (x,y)  ■  -y  +  x  +  1,  y<0)  -  1,  using  numdivx  -  10. 
dx 

TABLE  4.  Function  Library  Input  for  Diffel  Program, 

FUNCTION  LIBRARY  Model  Name:  diffel 

Name:  Expression: 

f  1-y+x; 

initl  initlx"0;endx-l;nuailntvl-10;lnltly"l; 


2.  Modified  Euler's  method: 

The  modified  Euler's  method  is  a  predictor-corrector  method. 

The  difference  equation  associated  with  the  method  is: 

yi  “  yi-l  yi-i>  +  f{*i»yi-l  +  “W-l^i-l^1  +  Error, 

2 

Error  -  0(h2). 

i  ■  l,2,.,,n,  where  n  »  number  of  intervals  and  x4  -  xQ  +  ih. 
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The  user  shall  enter  the  input  information  in  the  "function  library" 
of  Micro  SAINT  as  shown  in  Euler's  method* 

3 .  Runga-Kutta  method  of  order  four : 

Rvngui-Kutta  method  of  order  four  is  a  high  accuracy  method, 
Errors  -  O(h^),  but  it  requires  more  computation  pot  step.  The  difference 
equations  associated  with  the  method  are: 

kl  "  “<*1-1*  y1.l), 

k2  -  hf(xi_1  +  h/2,  yi_1  +  k/2), 

k3  -  hfU^  +  h/2,  yi-1  +  k/2) , 
k4  -  hf(xif  yA_ L  +  k3), 

yi  "  yi-l  +  <kl  +  2k2  +  2k3  +  k4^6» 

i  m  1,2,*** ,n • 

The  user  shall  enter  the  input  information  in  the  "function  library"  of 
Micro  SAINT  as  shown  in  the  example  in  Table  3* 

b.  First-order  system  of  ordinary  differential  equations: 


We  consider  two  differential  equations  in  two  unknowns  u^  and  u?: 

3xl  ' 

-~2  »  f2(x,ul,u2),  u1(xQ)  *(J  u2(xq)  m/3  • 


The  difference  equations  associated  with  the  extension  of  Runga* 
Kutta  method  to  systems  of  differential  equations  are: 


1,1  ■ 

hf 

i<xj 

,u1,u2) 

>  1  *  L|2| 

2,i  “ 

hf 

i(Xj 

+  h/2, 

Ujj  +  0»5lc^^fti2j  ^  ^ 

3,i  " 

hf 

i<*j 

+  h/2, 

t»ij  ^  0* ^^21>u2j  ^  ^  18 

M  “ 

hf 

i<xj 

+  h,  u 

l,j  +  k3i>u2j  +  k32^»  1  “  l*2» 

1,J+1 

- 

ul,j 

+  (ku 

+  2k2^  +  2k31  “*■  k4l)/^» 

2,j+l 

ai 

U2*-} 

+  (klf 

2  2k22  +  2k32  +  k42^^* 
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As  an  illustration  of  an  application,  we  use  an  example  [2],  about 
the  use  of.  Kirchkoff's  Law  in  circuit  theory.  Assuming  that  the  switch  in 
the  circuit  shown  in  Fig.  1  is  closed  at  time  t  ■  0, 


4  ohms 


2  henrys 


Figure  1.  RLC  electrical  circuit* 


the  currents  l^(t)  and  l2(t)  in  the  left  and  right  loops,  respectively,  are 
solutions  of  the  following  equations: 

|il  -  f1(t,I1,I2)  -  —4 1 jl  +  3I2  +  6, 1^(0)  -  0, 

~2  -  f2<t,ll,I2)  -  -2.4I1  +  1.6I2  +  3.6,I2(0)  -  0. 

The  user  shall  enter  the  input  information  using  x,  u1,  and  u2,  re¬ 
spectively,  for  t,  Ip  and  I2  in  the  "function  library"  of  Micro  SAINT  as 
shown  in  Table  5. 

TABLE  5.  Function  Library  Input  for  Diffe4  Program. 


FUNCTION  LIBRARY  Model  Name:  diffe4 

Name:  Expression: 

fl  (-4)*ul+3*u2+6; 

f 2  (-2.4)*ul+1.6*u2+3. 6; 

lnitl  initlx*'0;endx"l;numintvl-10;initlul-0;initlu2'-0 
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c.  Second-order  ordinary  differential  equation; 


To  approximate  the  solution  of  a  general  second-order  ordinary 

differential  equation 

d2y  -  f2(x,y,d£),  a  <  x  <  b, 
dx^  dx 

y(a)  -<2  ,  djr(a)  “j9, 
dx 

one  needs  to  transform  the  equation  into  a  system  of  first  order  using  the 

transformation  =  y,  u2  “  d_£>  to  gets 

dx 


iili  - 

dx 

du2  «. 

dx 


u2» 

^2^x,ul,u2^ * 


u1(a)  ,u2(a)  -  . 

To  illustrate  the  method,  we  consider  the  problem 

d2v  -  f (x,y ,dy)  -  2y'/x  -  2y/x2  +  xln(x) 
aJ*  A ^ 


dx 

y(i)  -  l, 

The  transformed  system  is* 

-  u2» 


dx 


y  (1)  -  0. 


dx 

du 

dx 


2  »  2u2/x  -  2u^/x2  +  x2ln(x), 


ut(l)  -  1,  u2(i)  -  0. 

The  user  ueeds  to  enter  the  initial  conditions  and  the  functional  form  of 
f2(x,ultu2)  as  shown  in  Table  6,  where  y  is  replaced  by  u^  and  dy/dx  is 
replaced  by  u2. 
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TABLE  6.  Function  Library  Input  for  Diffe5  Program. 

FUNCTION  LIBRARY  Model  Name:  diffe5 

Name:  Expression: 

f  2  2*u2/x+(-l)*2*ul/x''2+x*ln(x) 

ini tl  ini tlx-I ;endx-3 ;numin  tvl“40 { ini tiul-l ; ini tlu2-0 


4.  CONCLUSIONS 

We  developed  Micro  SAINT  programs  for  numerical  methods  of  integration 
and  differential  equations.  These  computational  modules  could  be  used  as 
subnetworks  in  modeling  problems  of  interest  in  naval  aerospace  medical 
research.  There  are  many  other  numerical  methods  for  Integra -ion  and  for 
solving  differential  equations.  Many  of  these  techniques  can  be  programmed 
in  Micro  SAINT*1  despite  its  minor  shortcomings,  such  as  the  absence  of 
defined  values  for  the  base  of  natural  logarithm,  and  the  capability  of 
reading  input  data  files. 

The  numerical  integration  methods  presented  here  are  adequate  when  the 
function  being  evaluated  is  relatively  simple,  that  is,  does  not  require 
many  time-consuming  manipulations.  The  differential  equations  methods, 
especially  the  Runga-Kutta,  are  adequate  for  problems  where  the  function  is 
easy  to  evaluate  and  the  accuracy  needed  is  small  (about  10"^  for  Runga- 
Kutta  methods). 

For  further  details  and  recommendations  on  which  method  to  use  for 
solving  a  given  nonstiff  initial-value  problem  (i.e.,  nonstiff  differential 
equations  with  initial  conditions),  we  recommend  that  the  papers  by  Hull  et 
al.  [3],  and  Enright  and  Hull  [4]  be  consulted.  For  details  on  methods  for 


*R,R.  Stanny,  Naval  Aerospace  Medical  Research  Laboratory,  Pensacola, 
FL,  personal  communication,  February  1989. 
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stiff  differential  equations,  we  recommend  consulting  Gear  {5),  Lambert 
[6],  Shampine  and  Gear  [7],  or  Enright  et  al.  i.8], 
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APPENDIX  A 


Micro  SAINT  Diagrams  and  Computer  Programs 


APPENDIX  A 


In  this  appendix,  we  give  actual  Micro  SAINT  diagrams  and  computer 
programs  for  the  following  algorithms: 

a.  Numerical  integration  of  a  function  f(x)  using  Trapezoidal  Rule. 

b.  Numerical  integration  of  a  function  f(x)  using  Simpson's  Rule. 

c.  Composite  Simpson's  Rule  for  double  integration  of  a  function 

f(x,y). 

d.  Composite  Simpson's  Rule  for  triple  integration  of  a  function 
f(x,y,z), 

"S 

e.  Euler  method  for  solving  a  first-order  ordinary  differential 
equation. 

f.  Modified  Euler  method  for  solving  a  first-order  ordinary  differen¬ 
tial  equation. 

g.  Runga-Kutta  method  of  order  four  for  solving  a  first-order 
ordinary  differential  equation. 

h.  Runga-Kutta  method  of  order  four  for  solving  a  first-order  system 
of  ordinary  differential  equations. 

i.  Runga-Kutta  method  of  order  four  for  solving  second-order  (linear 
or  nonlinear)  ordinary  differential  equations. 


A-l 


Figure  2.  Micro  SAINT  diagram  for  Trapezoidal  Rule  for  integrating  a  function  f(x). 


TASK  NETWORK  FOR  INTGTRPZ  PROGRAM 


0>  Types 


1  i 

start/trapez 
be  task  or  network): 

Type: 
Task 
Task 
Task 


(2)  Type: 


Network  Number:  0 
(1)  Name:  intgtrpz 

(3)  Upper  Network: 

(4)  Release  Condition: 

(5)  First  sub-job:  1 

(6)  Sub-jobs  (each  can 

Number :  Name : 

1  start/trapez 

S  addarea 

3  end 

Task  Number:  1 

(I)  Name:  start/trapez 

(3)  Upper  Network:  0  intgtrpz 

(4)  Release  Condition:  lj 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0$ 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect: 

initl ; j=l ;sum=0;delx=<endx-ini tlx) /numintvl ;x=initlx; 
f ;prevf=f ; 

(10)  Decision 
Following 

Number : 

(II)  a 

(13) 

(15) 

(17) 

(19) 

(El) 

(£3) 


Type:  Single  choice 
Task/Network : 

Name: 

addare  (la) 
(14) 
(16) 
(18) 
(SO) 
<BB) 
(34) 


Probability 
This  Path: 
15 


Of  Taking 


Task  Number:  2 
(1)  Name:  addarea 

(3)  Upper  Network:  0  intgtrpz 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0; 

(8)  Task's  beginning  effect: 

<9)  Task’s  ending  effect:  x=initlx+j#delx;f } 


(8)  Type: 


Task’s  ending  eTTect:  x=in 
sum-sum+(prevf+f ) ;prevf=f } j=j+lj 
(10)  Decision  Type:  Tactical 
Following  Task/Network ■ 


(11) 

(13) 

(15) 

(17) 

(19) 

<an 

(83) 

Task 

(1) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


(11) 
(13) 
(15) 
(17) 
(19) 
(21 ) 
(23) 


Number : 
2 
3 


Name: 

addare 

end 


Tactical  Expression: 


<=  numintvl; 
>  numintvl; 


Normal 


(18) 

(14) 

(16) 

(18) 

OO) 

(22) 

04) 


Number :  3 

Name:  end 

Upper  Network:  0  intgtrpz 
Release  Condition:  1; 

Time  Distribution  Type: 

Mean  Time:  0; 

Standard  deviation:  0; 

Task’s  beginning  effect: 

Task’s  ending  effect:  sum=sum*delx/2; 
Decision  Type:  Last  task 
Following  Task/Network; 

Number :  Name : 

( IE) 

(14) 

(16) 

(18) 

OO) 

(22) 

(24) 


(2)  Type: 


Probability 
This  Path: 


Of  Taking 


A- 3 


Network 


Task 


Task 


Task 


TASK  NETWORK  FOR  INTGSIMP 
Network  Number:  0 

(1)  Name:  intgsimp  (2)  Type:  Network 

(3)  Upper  Network: 

(4)  Release  Condition:  1; 

(5)  First  sub-job:  1  start/simpson 

<65  Sub- jobs  (each  can  be  task  or  network): 

Number:  Name?  Type: 

1  start/simpson  Task 

2  adoorea  Task 

3  end  Task 


Task  Number:  1 

(1)  Name:  siart/simpscn 

(3)  Upper  Network:  0  intgsimp 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

<75  Standard  deviation:  0} 

(8)  Task's  beginning  effect: 

(9)  Task’s  ending  effect:  initlj 
delxMendx-initlxf/numintvl  im^numintvl/Sj 
j=l ; sum=0  »xin=initlx!x=xin»fi tempf “f » 

(10)  decision  Type:  Single  choice 

Following  Task/Network:  Probabi 

Number:  Name:  This  P 

<115  2  addare  (12)  1; 


Type: 


Probability  Of  Taking 
This  Path: 
li 


Task  Number: 


(I)  Name:  addarea  (2)  Type:  Task 

<3)  Upper  Network:  0  intgsimp 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 
<6)  Mean  Time:  0; 

(7)  Standard  deviation:  0} 

(8)  Task's  beginning  effect: 

(9)  Task's  ending  effect:  x=xin+delx;f }f l=f »x=x+delx Jf ;f2=f ; 
sum*sum+tempf +4#f  f +f 2 ; 
xin=xj 

xin*x;tempfaf2; j=j+l ; 

(10)  Decision  Type:  Tactical 

Following  ■ ask/Network :  Tactical  Expression: 

Number:  Name: 

(II)  2  addare  (12)  j  <  m+1; 

<  13)  3  end  (14)  j  >=  m+1! 

(15)  (16) 

(17)  (18) 

(19)  (20) 

(21)  (22) 

(23)  (24) 


Tactical  Expression: 


Task  Number:  3 

(I)  Name:  end  (2)  Type:  Task 

(3)  Upper  Network:  0  intgsimp 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0{ 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect:  sum=sum*delx/3i 

(10)  Decision  Type:  Last  task 

Following  Task/Netwdrk:  Probability  Of  Taking 

Number:  Name:  This  Path: 

(II)  ( 12) 


A- 5 


TASK  NETWORK  FOR  DBLINTEG  PROGRAM 


Network  Number;  A 
(1)  Nam?:  dblinteg 

(3)  Upper  Network: 

(4)  Release  Condition:  1; 

<5)  First  sub-jobi  1  start 

(6)  Sub-jobs  (eacU  can  be  task  or  network): 
Number:  Name:  Typo: 

1  start  Task 

8  integxl  Task 

3  integy  Task 

4  integxS  Task 

5  end  Task 


(8)  Type:  Network 


Task 
(1) 

<3) 

(4) 

(5) 

(6) 

(7) 

<C) 

(9)  Task's  ending 
turn 1=0} 5umS=0 } sum3=0 i i =0 j 


dblinteg 

1! 

Type:  Normal 


(8)  Type:  Task 


Number :  1 

Name:  start 
Upper  Network:  0 
Release  Condition: 

Time  Distribution 
Mean  Time:  0; 

Standard  deviation:  0; 

Task’s  beginning  effect: 

Task’s  ending  effect:  initl }delx=(endx-initlx) /numdivxj 


(10) 


(11) 
(13) 
(15) 
(17) 
<  19) 
(81) 
(83) 


Decision  Type:  Single  choice 
Following 
Number : 

8 


Task /Net work: 
Name : 
integx 


(18) 

(145 

(16) 

(IB) 

(80) 

(88) 

(84) 


Probability 
This  Path; 
i; 


Of  Taking 


Task  Number:  8 

(1)  Name:  integxl  (8)  Type:  Task 

(3)  Upper  Network:  0  dblinteg 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0} 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect:  x=initlx+i*delx;funcdofx;funccofx; 
dofx=funcdofxicofx-funccofx;HX=(dofx-cofx)/numdivy; 
y=dofx;f ;f l=f ;y=cofx;f ;f3=f jtempi«8#int( i/E) ; 

kl=f l+fa;kB-0;k3=0; 


(11) 

(13) 

(15) 

(17) 

(19) 

(81) 

(S3) 


Decision  T 
Following 
Number : 

o 


ype:  Single  choice 
Task/Network: 

Name: 

integy  (IS) 
(14) 
(16) 
(18) 
(SO) 
(88) 
(84) 


Probability  Of  Taki 
This  Path: 
i; 


ng 


Task  Number:  3 

(I)  Name:  integy  (8)  Type:  Task 

(3)  Upper  Network:  0  dblinteg 

(4)  Release  Condition:  li 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  Oj 

(7)  Standard  deviation:  0; 

(B)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect:  y=cofr+j#HX;f ;z=f ; 
temp=8*int( j/S) ; 

if  temp  ==  j  then  kE=kP^z  else  k3=k3+z> 

L=(H+8*kS+4*k3)#HX/3; 

tur  Decision  Type:  Tactical 

Following  Task/Network:  Tactical  Expression: 

Number:  i 

(II)  3  integy  (18)  j<*numd ivy-1 S 

(13)  4  integx  (14)  j  >  numdivy-lj 


A- 7 


Task  Number: 


Name:  integxE  (8)  Type:  Tast 

Upper  Network:  0  dblinteg 
Release  Condition:  i; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0; 

Task's  beginning  effect: 

Task's  ending  effect: 

==  0  I  i  «■  numdivx  then  suml*suml+L  else  if  tempi==i  then 


sum2=sumE+L  else  sum3=sum3+L 5 
if  i  ==  numdivx  then  endinteb=l! 
(10>  Decision  Type:  Tactical 
Following  Task/Network: 

kl L  .  .  .  k.  I 


Number : ' 

a 

5 


Name : 

integx 

end 


i=i  +  l  * 

Tactical  Expression: 

endinteg==0; 
endinteg==l ; 


Number :  5 

Name:  end  (E)  Type:  l 

Upper  Network:  0  dblinteg 
Release  Condition:  l; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0; 

Task’s  beginning  effect: 

Task’s  ending  effect:  sum=<suml+B#sum2+4#sum3)*delx/3; 
Decision  Type:  Last  task 


Decision  Type:  Last  ta< 
Following  Task/Network? 
Number:  Name: 


Probability  Of  Taking 
This  Path: 
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TEST  NETWORK  FOR  TRIPLINT  PROGRAM 


Network  Number:  0 
(1)  Name:  triplint 

(3)  Upper  Network: 

(4)  Release  Condition: 

(5)  First  sub-job;  1 

(6)  Sub-jobs  veacn  can 

Number :  Name : 

1  start 

2  integl 

3  integ3 

4  integ4 

5  integ5 

6  integ6 

7  end 


Type:  Network 


start 
be  task 


network) 
Type : 
Talk 
Task 
Task 
Task 
Task 
Task 
.  Task 


!ask 

Jj 

I 

(6) 

(7) 

(8 

(9 

ini  tl 


Number :  1 

Name:  start  (2)  Type:  Task 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0: 

Task  s  beginning  effect: 

Task  s  ending  effect: 

:  n=numdivx/2 : m=numdivy/2 ; p-numdivz/2 ; h= (endx- ini  tlx) / (2#n) ; 
j2  =  0;j3  =  0:i  =  <&;j=0;  *  ’  . 

Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 

2  integl  (12)  1 ; 


k  1  =  0 ; 
(10) 


Number :  2  „  v  „ 

Name:  integl  (2)  Type:  Task 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0; 

Task.s  beginning  effect:  .  .  .  ,  . 

Task’s  ending  effect:  x« mi t lx+ i #h ; dx ; cx ; hx= (dx-cx) / ( 2#m) ; 


Mean  Time:  v; 
Standard  deviation: 
Task.s  beginning  effe 
Task  s  ending  effect: 


k2  =  0;k3  =  0: j  =0 ; 

Decision  Type:  Single 
Following  Task/Network; 
Number :  Name ; 

3  integ3 


cho ice 


Probability  Of 
This  Path: 

1 ; 


Taking 


Number:  3 

Name:  integ3  12)  Ty 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0: 

Task.s  beginning  effect: 

Task  s  ending  effect: 

■  j#hx : betaxy ; zfbetaxy ; f xyz ; f l  =  f xyz ; alnhaxy ; 
jetaxy-alphaxy) / (2#p; ; z=alphaxy : f xyz ; J2=f xyz ; Ll=f 1+ 
;  L3  =  0 ; ks 1 ; 

Decision  Type:  Single  choice 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  v  This  Path: 

4  integ4  (12)  1  ; 

14) 

(16) 


Type:  Task 


f  l  +  f2; 


A- ’10 


4 

5 

(6) 

7) 

8 

(9) 
if  i 
k=k  + 

(10) 


Number  :  ,  4 

Name:  integ4  (2)  Typ 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time :  0 ; 

Standard  deviation:  0; 

Task;s  beginning  effect: 

Task's  ending  effect:  z=alphaxy+k#hy  j  f  xyz  ;  qq  =  f  xyz  ; 
nt  (k/2)  ==k/2  then  L2  =  I.2  +  qq  else  L3=L3+qq; 

decision  Tvoe :  Tactical 

Following  Task/Network:  Tactical  Expression:.' 


Type:  Task 


Number : 

5 

4 


Name 

intej 

intej 


Tactical  Expression: 

k  =  =  2  #p ; 
k  <  2#p; 


Task 
( 1) 

3 

4 

5 

6 
7 
0 

k3  =  ^3 

\itil 


Number:  5 

Name:  integ5 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0: 

Task  s  beginning  effect: 

Task  s  ending  effect:  L=(Ll+2# 
=0  !j==2  then  kl=kl+L  else  if  i 

i  +  L ; 

decision  Type:  Taatical 
Following  Task/Network:  Ta 

Number :  Name : 

6  integ6  12)  j= 

3  inte|3  (14)  j 


Type:  Task 


lseLif 2intt j /2?== j^23then  k2=k2+L  else 


Tactical  Expression 

j==2#m+l ; 

3  <2#m+l; 


Task 
(  1) 

3 

(4 

5 

6 

7 

8 
(9 

if  i 
else 


Number :  6 

Name:  integ6  (2)  Type:  Task 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0: 

TaSk.S  beginning  effect:  ,  ,  . 

Task  s  ending  effect:  ksum- (kl +  2#k2t4#k3) #hx/3 :  .  ■  _  . 

=  =  0  !  i*=2#n  then  jl=jl+ksum  else  if  int (i/2) =-i/2  then  j2*j2+ksum 
j3=j3+ksum; 


decision  Typ 
Following  Ta 
Number : 

7 

2 


pe :  Taatical 

ask/Network : 
Name : 
end  ( 

integl 


Tactical  Expression: 

i==2#n+ 1 : 
i  <2#n+l; 


Task 
$  1) 

3 

4 

5 

6 

8 

9) 

( 10) 


Number:  7 

Name :  end 

Upper  Network:  0  triplint 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0; 

Standard  deviation:  0: 

Task  S  beginning  ef feet : 

Task  s  ending  effect:  jsura-(jl' 


Type:  Task 


Decision  Type:  Last  ta 
Following  Task/Network: 
Number:  Name: 


feet:  j sum- ( j 1+ j 2*2+4# j 3) #h/3 ; 

Last  task 

etwork:  Probability  Of  Taking 

ame :  This  Path: 


A-ll 


TASK  NETWORK  FOR  DIFFE1  PROGRAM 
Network  Number:  0 

(1)  Name:  diffel  (2)  Type:  Network 

(3)  Upper  Network: 

(4)  Release  Condition:  1; 

(5)  First  sub-job:  1  start/  Euler 

(6)  Sub- jobs  <each  can  be  task  or  network): 

Number:  Name:  Type: 

1  start/  Euler  Task 

2  Euler  Task 

3  end  Task 


Task  Number:  1 

(1)  Name:  start/  Euler  (2)  Type:  7 as 

(3)  Upper  Network:  0  diffel 

(4)  Release  Condition:  1} 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0; 

(S)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect:  initl;delx=(endx-initlx)/numintvl ; 
xin=initlxjyout=initly ; 
k=>l  5 

< 10)  Decision  Type:  Single  choice 


Following 
Number : 

2 


ype:  Single  c 
Task/Network : 
Name: 
Euler 


Probability  Of  Taking 
This  Path: 

15 


Task  Number:  2 

U)  Name:  Euler  (2)  Type:  Task 

(3)  Upper  Network:  0  diffel 

(4)  Release  Condition:  1} 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0; 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect: 

x=x  in  5  y=yout  ;f  5yout=yout+del  x#f  *  k-k+1  *  x  in=x in+del x ; 

(10)  Decision  Type:  Tactical 


Following 

Task/Network: 

Tactical  Expression: 

Number : 

Name: 

(11) 

2 

Euler  (12) 

kOnumintvl » 

(13) 

3 

end  (!4) 

k>numintvl 5 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

Task  Number:  3 

(1)  Name:  end  (2)  Type:  Task 

(3)  Upper  Network:  0  diffel 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  1 5 

(7)  Standard  deviation:  0; 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect: 

(10)  Decision  Type:  Last  task 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 


(11) 


(12) 
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TASK  NETWORK  FOR  DIFFE2  PROGRAM 


Network  Numbers  0 
(1)  Names  diffeE 

(3)  Upper  Network: 

(4)  Release  Conditii 


Types  Network 


First  sub-jobs 
Sub-jobs  (each 


Number  s 
i 

a 

3 

Task  Numbers 


-jobs  1  start/mod.  Euler 
(each  can  be  task  or  network): 
Name:  Type: 

start/mod.  Euler  Task 


modif-Euler 

end 


Task 

Task 


Name:  start/mod.  Euler  (2)  Type:  Ta* 

Upper  Networks  0  diffeE 
Release  Condition:  1; 

Time  Distribution  Type:  Normal 
Mean  Time:  0  j 
Standard  deviation:  0; 

Task's  beginning  effect: 

Task's  ending  effect:  initl>delx=(endx-initlx )/numintvl ; 


(B)  Task's  beginning  effect: 

(9)  Task's  ending  effect:  ini t 1 » 
xinsinitlxjyout^initly; 

k=i; 

(10)  Decision  Type:  Single  choice 
Following  Task/Network: 

Numbers  Name: 

(11)  a  modif-  (IB) 


Probability  Of  Taking 
This  Paths 
i; 


Task  Number:  E 

(I)  Name:  modif-Euler  (E)  Type:  Task 

(3)  Upper  Network:  0  diffeE 

(4)  Release  Conditions  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0} 

(8)  Task's  beginning  effect: 

(9)  Task’s  ending  effect:  x=xin}y«yout;f {y=yout+delx*f }f 5cf=f ? 
yout=yout+delx*(f+cf )/E;k=k+l } 

xin=xin+delx! 

(10)  Decision  Type:  Tactical 

Following  Task/Network:  Tactical  Expression: 

Number:  Name: 

(II)  E  modif-  (IE)  k<-numintvl» 

(13)  3  end  (14)  k>numintvl; 

(15)  (16) 

(17)  (18) 

(19)  (SO) 

(81)  (EE) 

(23)  (E4) 


Task  Number:  3 

(1)  Name:  end  (E)  Type:  Task 

(3)  Upper  Networks  0  diffeE 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  1? 

(7)  Standard  deviation:  0! 

(8)  Task's  beginning  effect: 

(9)  Task’s  ending  effect: 

(10)  Decision  Type:  Last  task 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 

(12) 


(11) 
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Figure  8.  Hicro  SAINT  diagram  for  Runga  Kutta  method  for  solving  a  differential  equation. 


TASK  NETWORK  FOR  DIFFE3  PROGRAM 


Network  Numbers 


Names  diffe3 
Upper  Networks 
Release  Conditions 
First  sub- job s  1 
Sub- jobs  (each  can 


(2)  Types  Network 


Number  s 
1 
2 
3 

Task  Numbers 


Conditions  1; 
a- job s  1  start/  R-K-4 
(each  can  be  task  or  network )s 
Names  Types 

start/  R-K-4  Task 

runga-kutta  4  Task 

end  Task 


Names  start/  R-K-4  (2)  Types  Tat 

Upper  Networks  0  diffe3 
Release  Conditions  1; 

Time  Distribution  Types  Normal 
Mean  Time:  0; 

Standard  deviations  0; 

Task’s  beginning  effects 

Task’s  ending  effects  initl ?delx=(endx-initlx)/numintvl ! 


xin=initlx}yout=initly;ksl » 

(10)  Decision  Types  Single  choi 
Following  Task/Networks 


Allowing 
Number  s 
2 


Names 

runga- 


Probability  Of  Taking 
This  Paths 
i; 


Task  Numbers  2 

(1)  Names  runga-kutta  4  <2 

(3)  Upper  Networks  0  diffeS 

(4)  Release  Conditions  1; 

(5)  Time  Distribution  Types  Normal 

(6)  Mean  Times  0} 

(7)  Standai  d  deviations  0; 

(8)  Task's  beginning  effects 

(9)  Task’s  ending  effect: 

x=x in ; y^yout ; f ; k l=del x#f ; x=x in+de J  x/2 }y=yout+k 1 /2  s 
f ;  k2=de  1  x#f ;  y^yout+k2/2 }  f }  k3sde  1  r,#f ;  x»x  in+del  x ; 
y=yout+k3;f 5k4=delx#f ; 
yout=yout+(kl+2*(k2+k3)+k4)/6i 
xin=xjk“k+l; 


Types 


Decision  Typ 
Following  Ta< 
Number : 

2 

3 


pes  Tactical 
ask/Networks 
Names 

runga-  <  I 
end  ( 1 


Tactical  Expressions 

k<*r.umintvl ; 
k>numintvl 5 


Task  Numbers  3 

(I)  Names  end  (2)  Types  Task 

(3)  Upper  Networks  0  diffe3 

(4)  Release  Conditions  1; 

(5)  Time  Distribution  Types  Normal 

(6)  Mean  Times  l; 

(7)  Standard  deviations  0; 

(8)  Task's  beginning  effects 

(9)  Task's  ending  effects 

(10)  Drcision  Types  Last  task 

Following  Task/Networks  Probability  Of  Taking 

Numbers  Names  This  Paths 

(II)  (12) 
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TASK  NETWORK  FOR  DIFFE4  PROGRAM 


Network  Number: 

<1)  Name:  diffe4 
(3)  Upper  Network: 

<4)  Release  Condition: 
<5)  First  sub-job:  1 
<6)  Sub- jobs  <each  can 
Number :  Ni 

1  si 

2  s' 

3  er 

Task  Number:  1 


(2)  Type:  Network 


Condition:  1 ; 

3- job:  1  start/  system 
<each  can  be  task  or  network): 
Name:  Type: 

start/  system  Task 

systm/lin/nonlin  Task 

end  Task 


Name:  start/  system  <2)  Type:  Ta« 

Upper  Network:  0  diffe4 
Release  Condition:  l; 

Time  Distribution  Type:  Normal 
Mean  Time:  Oj 
Standard  deviation:  0? 

Task's  beginning  effect: 

Task’s  ending  effect:  initi;delx3<endx-initlx)/numintvl; 
ini tlx?ulout=initlul  ;u2out=initlu2; 


^  /  t  UWOllUai  UCYiOViU  II*  v  j 

<0)  Task's  beginning  effect: 

(9)  Task’s  ending  effect:  ini t ] 
xin=ini tlx>ulout=initlul  ;u2out=i 


Decision  Type:  Single  choice 
Following  Task/Network: 


Following 
Number : 

2 


fask/Network : 
Name: 
systm/ 


Probability  Of  Taking 
This  Path: 
i; 


Task  Number: 
< 1 )  Name : 
<3)  Upper  N 


Numoer ;  e 

Name:  systm/lin/nonlin 
Upper  Network:  0  diffe4 
Release  Condition:  1; 

Time  Distribution  Type:  t* 
Mean  Time:  0; 

Standard  deviation:  0; 
Task’s  beginning  effect: 
Task’s  ending  effect:  x=> 


Type: 


Normal 


a 

x=xin!ul=ulout ;u2=u2out i 


uiout=uiout+<kH+2*:k2i+kai  )/o; 
u2out=u2out+ ( k 12+2# ( k22+k32 )  +k42 ) /6 ! 
xinsxin+delxik^k+l J 


ul.uu  v  -uuuu  n  \  rs 

xinsxin+delxik^k+l J 
(10)  Decision  Type:  Tact 
Following  Task/Netwo 

w  i,  .  ^ : ki . 


Following 
Number : 

2 

3 


pe:  Tactical 
ask/Network : 
Name: 
systm/  ( 
end  < : 


Tactical  Expression: 

k<=numintvl * 
k>numintvl J 


Task  Number:  3 

<1)  Name:  end  <2)  Type:  Task 

<3)  Upper  Network:  0  diffe4 
<4)  Release  Condition:  1; 

<5)  Time  Distribution  Type:  Normal 
<6)  Mean  Time:  1} 

(7)  Standard  deviation:  0; 

<0)  Task’s  beginning  effect: 

<9)  Task’s  ending  effect: 

UO)  Decision  Type:  Last  task 

Following  Task/Network:  Probability  Of  Taking 

Number:  Name:  This  Path: 

<11)  (12) 
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TASK  NETWORK  FOR  DIFFE5  PROGRAM 


Network  Numbers  0 
(1)  Names  diffe5 
<3)  Upper  Networks 

(4)  Release  Condition: 

(5)  First  sub-job?  1 


<£)  Type:  Network 


i; 

start/End 


order 


Numbers 

1 

a 

3 


Sub- jobs  (each  can  be  task  or  network): 


Name: 
start/End 
End  order 
end 


order 

1 in/nonl in 


Type: 

Task 

Task 

Task 


Task  Number: 

(1)  Name:  start/End  order  <£)  Type:  Tas 

(3)  Upper  Network:  0  diffeS 

(4)  Release  Condition:  1; 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0; 

(E)  Tisk's  beginning  effect: 

(9)  Task's  ending  effect:  initl jdelx=(endx-initlx)/numintvl > 
xin^initlxjulout^mitlul  ;u2out=initlu2; 

k=l ;error=0; 

(10)  Decision  Type:  Single  choice 


Type: 


Following 
Number : 

a 


ype:  Single  c 
Task/Network: 
Name: 

End  or  1 


Probability  Of  Taking 
This  Path: 
i; 


Task  Number: 


(1)  Name:  End  order  1 in/nonl in  (E)  Type:  Ti 

(3)  Upper  Network:  0  diffe5 

(4)  Release  Condition:  1? 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  Time:  0; 

(7)  Standard  deviation:  0> 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect:  x=xin;ul=ulout;u£=u£out jf l=uEjfE! 
k  1 1 =de 1 x*f 1 ; k 1 E=del x*f 2 ; x-x i n+de 1 x /E » 
ul=ulout+kll/2;u2=u£out+ki£/25fl=uE5fe; 
kSl=delx#fi;k£2=delx#f2;ul=ulout+k21/£; 
u2=u£out+k2E/2;f I=u£;f2;k31=delx*f 1 5k3E=delx#f£! 
ul=ulout+k3i ;uE=uEout+k325x=xin+delx»f l=uEjfE{ 
k41=delx*fl;k4E=delx*f2i 
ulout=ulout+(kll+£*(k£l+k31 )+k41)/65 

u£out-uScu  t+  <  k  1E+E*  ( kSE+k3E )  +k4£ ) /6i 
x in-xin+del x ;k=k+l * 

(10)  Decision  Type:  Tactical 


Following 

Task/Network: 

Tactical  Expression 

Number : 

Name: 

(11) 

2 

End  or  (IE) 

k<=numintvl * 

<  13) 

3 

end  ( 14 ) 

k>numintvl ; 

(15) 

\  16) 

(  17) 

(IB) 

(  19) 

(SO) 

(£1) 

(EE) 

(S3) 

(E4 ) 

Task  Number:  3 

(I)  Name:  end 

(3)  Upnor  Network:  0  diffe5 

(4)  Release  Condition:  1* 

(5)  Time  Distribution  Type:  Normal 

(6)  Mean  \ ime:  0! 

(7)  Standard  deviation:  0* 

(8)  Task’s  beginning  effect: 

(9)  Task’s  ending  effect: 

(10)  Decision  Type:  Last  task 

Following  Task/Network:  Pr 

Ni  rnber :  Name :  T 

(II)  <1S) 


Type: 


Probability  Of  Taking 
This  Path: 
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APPENDIX  B 


Graphical  Representation  for  Solutions  of  Differential  Equations. 


APPENDIX  B 


In  this  appendix  we  present  actual  Micro  SAINT  graphical  outputs  for 
solutions  of  the  following  cases: 
a*  Runga-Kutta  method  for 

y'  at  1-y  +  x>  y(0)  -  1. 

b.  Runga-Kutta  method  for  the  system: 

£IX  -  f1(t,I1,I2)  -  -4Il  +  3I2  +  6,^(0)  -  0, 
dt 

dl2  “  f2^*Il*I2^  "  “2*4I1  +  1«6I2  +  3.6,I2(0)  -  0. 
dt 

c.  Runga-Kutta  method  for  second-order  differential  equation: 
x2y"  -  2xy"  +  2y  -  x3Ln(x),y(l)  -  l,y'(l)  -  0, 


B-l 


